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An increasing number of measurements in fundamental and applied physics rely on magnetically shielded 
environments with sub nano-Tesla residual magnetic fields. State of the art magnetically shielded rooms 
(MSRs) consist of up to seven layers of high permeability materials in combination with highly conductive 
shields. Proper magnetic equilibration is crucial to obtain such low magnetic fields with small gradients in 
any MSR. Here we report on a scheme to magnetically equilibrate MSRs with a 10 times reduced duration of 
the magnetic equilibration sequence and a significantly lower magnetic field with improved homogeneity. For 
the search of the neutron’s electric dipole moment, our finding corresponds to a linear improvement in the 
systematic reach and a 40% improvement of the statistical reach of the measurement. However, this versatile 
procedure can improve the performance of any MSR for any application. 


I. INTRODUCTION 

Magnetically shielded rooms (MSRs) are environments 
where external electromagnetic field distortions at low 
frequencies are strongly damped and the earth or am¬ 
bient magnetic fields are strongly suppressed. MSRs 
are typically used in biomagnetioi, medical and funda¬ 
mental physics applications^. For these, both a strong 
suppression of time varying external fields of frequen¬ 
cies from mHz to kHz, as well as the resulting residual 
quasi static magnetic field (< 1 nT) and residual mag¬ 
netic field gradients (< 1 nT/m) are important. Shield¬ 
ing in the low frequency range of up to tens of Hz is 
usually achieved by using large amounts of highly per¬ 
meable alloys (e.g. Permalloy). However, already at 
the 1 Hz regime the conductivity of the shielding mate¬ 
rial plays a significant role due to induced currents over 
large areas of the surface. For shielding of radio fre¬ 
quency (RF) disturbances, an additional layer of highly 
conducting material is added. The reference facility for 
magnetic shielding is the Berlin Magnetically Shielded 
Room 2 (BMSR-2) at the Physikalisch-Technische Bun- 
desanstalt (PTB) Berlin^. With 7 layers of Mumetall^ 
and an additional aluminum layer for RF shielding it has 
a passive shielding factor (SF) of 75000 at 0.01 Hz. The 
new shield for the next-generation neutron electric dipole 
moment (nEDM) experiment at the Technische Univer¬ 
sitat Munchen (TUM) reaches even higher SFs exceed¬ 
ing 1000000 at 1 mHz, without any active measures to 
compensate for low-frequency drifts^. Materials used for 
shielding have a high permeability and a non-zero rema- 
nence and therefore can be easily magnetized by arbi¬ 
trary fields. Such fields are typically generated by (i) 
static external sources, (ii) internally applied fields for 
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e.g. in low-field NMR applications, (iii) magnetic fields 
emerging from adjacent shielding layers and (iv) mag¬ 
netic distortions, e.g. from construction materials of the 
shield itself. To achieve lowest possible residual fields 
and field gradients, the material has to be ’degaussed’ 
repeatedly. ’Degaussing’ is the commonly used term for 
the reduction of the remanent magnetization in magne¬ 
tizable materials. It is usually performed by applying 
a sinusoidal current with decreasing amplitude to coils 
wound around the shielding material inducing magnetic 
flux in a closed loop of shielding material. The starting 
amplitude has to be large enough to reach saturation of 
the magnetization everywhere within the material. At 
the end of the cycle, the amlitude has to reach zero as 
exactly as possible. Any DC offset or noise can worsen 
the quality of the degaussing procedure. Degaussing of 
MSRs has been developed to an unprecedented qualityi^, 
with achieved residual fields of < 500 pT and gradients 
of < 0.2 nT/m over the center cubic meter volume in¬ 
side walk-in MSRs. These small fields can be generated 
in a highly reproducible manner to within 10’s of pT, 
independent of the initial conditions. 

The relationship between residual field and shielding 
has also been well investigated^. The shielding factor of 
an MSR is given as the frequency dependent attenuation 
of external disturbances 
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A a denotes the field magnitude at a position measured 
before the presence of the MSR, Ai the magnitude at the 
same position with the MSR, u is the frequency of the ex¬ 
ternal disturbance and 0 is a phase shift between the ex¬ 
citation field and the measured field inside the MSR. As¬ 
suming a sufficient amount of shielding material, in the 
static case (uj = 0) the field within the MSR is dominated 
rather by the residual field from the shielding material, 
and not by the effect of external fields. Ultra-low resid¬ 
ual fields are achieved by effectively offering the magnetic 
field a path of minimal magnetic resistance around the 
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volume to be shielded. By minimizing the potential en¬ 
ergy of the magnetic domains within the shielding ma¬ 
terial in their surrounding field an equilibrium state is 
formed, which is stable in time for a static environment. 
Therefore degaussing does not per se remove the mag¬ 
netic field as the name suggests, but rather optimizes 
the magnetic field in the material with respect to the ex¬ 
ternal condition. For these reasons, we prefer the usage 
of ’magnetic equilibration’ instead of degaussing to refer 
to this process (naturally, spots of strongly magnetized 
material still have to be degaussed to reduce their mag¬ 
netization). One of the procedures to minimize resid¬ 
ual magnetic fields inside MSRs is described in Ref. 0- 
Here, the coil setup for equilibration consists of 12 coils 
for each layer of the MSR with one along each edge of 
the cube. Four of these coils are connected in series to 
form a closed loop around one spatial direction (see also 
Fig [lja) . The three spatial directions are then degaussed 
sequentially. A sinusoidal current with a peak ampli¬ 
tude of 21x5 Amperexturns and a frequency of 10 Hz is 
applied. The amplitude decreases over 2000 periods (‘cy¬ 
cles’) to a value as close to zero as technically achievable. 
Equilibration therefore takes 600 s for each layer. With 
this procedure residual fields below 1.5 nT and gradients 
smaller than 2 nT/m have been achieved. In Ref. 0 one 
of the main criteria for a reproducible procedure is the 
step size of the amplitude decrease. The difference in 
amplitude between two subsequent maxima of the sinu¬ 
soidal current should be as small as possible as this delta 
corresponds to an error in the residual field. Therefore it 
is suggested to use several thousand cycles. 

In this work we describe a new procedure for mag¬ 
netic equilibration, which is capable of producing ex¬ 
tremely small residual fields, comparable to the best ever 
achieved. However, the result is achieved within a much 
shorter time constant of only 30-50 s per layer, with sim¬ 
pler wiring of equilibration coils and with simpler equip¬ 
ment. In the previous procedure with a duration of about 
900 s, residual fields of < 0.7 nT and field gradients of 
0.3 nT/m were achieved in the MSR at TU MimcheiA 
The new configuration can produce the same or better 
values in only 150 s. Here, both layers of the MSR are 
equilibrated in a defined sequence. For MSRs with only 
one layer of magnetizable material, such a sequence can 
even be shorter, on the order of only 10 s with a resulting 
residual field of <10 nT over a cubic meter inside, making 
single layer shielding a highly competitive technical solu¬ 
tion at low cost. Our finding has consequences for several 
applications, including in fundamental physics where the 
duty cycle of experiments is affected by frequent mag¬ 
netic equilibration sequences (e.g. required due to fre¬ 
quent reversal of applied NMR fields inside). Another 
example would be applications in biomagnetism, where 
the door to an MSR may be opened frequently. 



FIG. 1. The setup of the equilibration coils for the old con¬ 
figuration, shown in a), and the new configuration in b). In 
a) the 4 coils to create a magnetic flux around the Z direction 
and their connections are shown. The new configuration in b) 
has the old Y and Z direction connected in a L shape so that 
a magnetic flux is created in both directions simultaneously. 



FIG. 2. Allan deviation^ calculated for data recorded with 
the fluxgate inside the MSR for about 20000 s. The inset 
shows the corresponding time data, shifted to 0 at t = 0 
on the left axis. Bx and Bz have then been offset by 0.4 nT 
and -0.4 nT respectively for clarity. The smooth change of 
the magnetic field over time correlates with a smooth, mono¬ 
tonic drop of 0.2° C in the environmental temperature. To 
attribute the correlation of field and temperature to a change 
of permeability is however more complex. 
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II. APPARATUS 

The MSR used for this measurement and also the pre¬ 
vious configuration of equilibration coils are described in 
Ref. @ (see also Fig |TJa) . The two layers of the MSR are 
equilibrated individually, first the inner layer, then the 
outer one and after that the inner one again. For each 
layer, we start with the X direction followed by the Z 
direction and then the Y direction (see also FigQ]). It 
was found that equilibration of one direction is not suffi¬ 
cient: subsequent equilibrations in the second and third 
direction each improved the result. Note the door is in 
the XZ-plane at negative y coordinate. Each equilibra¬ 
tion sequence starts with an amplitude that saturates the 
magnetizable material, in the case of the inner layer that 
is 9 x 7 Ampere turns, and then linearly decreases over 
1000 cycles to below noise level. A frequency of 10 Hz 
is used. The waveform for equilibration is generated in 
software on a PC, converted to an analog voltage by a 
16 bit digital-analog converter (DAC) with a sampling 
frequency of 20 kHz. This voltage is passed through an 
analog voltage divider to use the full range of the DAC 
and a 100 Hz lowpass filter to smooth out DAC steps. 
It is converted to a current in a power amplifier— and 
passed through a transformer to remove any DC offset 
before being distributed to the respective coil of the MSR 
by an automated relay switch. In this way magnetic equi¬ 
libration can be done completely automated, in a process 
that takes about 15 min. 

Field maps have been measured using a three-axis flux- 
gate sensor—. The sensor heads at the end of flying leads 
are mounted in a nonmagnetic structure so that they 
form a right handed coordinate system. Their signal is 
read out with a 16-bit analog digital converter (ADC). 
Each data point in the maps is an average over a measur¬ 
ing time of 1 s, where data has been recorded with a sam¬ 
pling frequency of 1 kHz. Several measurements taken at 
the same point deviate less than 0.1 nT from each other. 
To correct for drifts of fluxgate offsets, a calibration is 
done before, during and after each field map. The sensor 
is calibrated in the center of the room where the fields 
are smallest, to keep the influence of alignment errors as 
small as possible. To determine the offset of a sensor the 
field is measured at a position and then again at the same 
position with the sensor rotated by 180° to measure in 
the opposite direction. Assuming a linear drift between 
calibrations the offsets drift rates were < 10 pT/min. All 
maps were corrected for these offset drifts. 

Figure [5] shows the performance of the fluxgate mag¬ 
netometer using an Allan variance plot. To obtain the 
optimal performance, an integration time of about 100 s 
would be required. For shorter times, noise will have 
an increasing effect on the measured signals. For longer 
times the measurement will be dominated by drifts that 
could come from changing offsets or drifts of the read¬ 
out electronics. To keep the measurement time per point 
reasonably short, an integration time of 1 s was chosen. 
Here the Allan deviation is already below 30 pT. With 
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FIG. 3. Maps of the magnitude of the magnetic flux den¬ 
sity inside the MSR after magnetic equilibration in three XY 
planes at different heights. The points are 10 cm apart in 
each direction. The columns show: the previous coil config¬ 
uration (a), and the L shaped configuration ( (b) and (c) ), 
where the later were taken with several days between the two 
measurements. 


this and repeated offset calibration an overall accuracy 
of ~ 0.3 nT can be reached, so that low-noise fluxgate 
sensors are well suitable for this kind of measurements. 

After magnetic equilibration in the earth field without 
active stabilization, field maps were recorded for the re¬ 
spective configuration. Each map consists of 108 points 
on a 6 by 6 grid in 3 different XY planes (for coordinate 
system see Fig. [T|) in the center of the MSR, 35 cm be¬ 
low and 35 cm above it. The distance between the grid 
points along each axis is 10 cm. The sensor was moved 
by hand from one point to the next with an estimated 
accuracy of about 0.5 cm and a tilt-precision of better 
than few degrees. The distance of the individual sensors 
in the mount is 2 cm, so the overall position accuracy of 
the measured points is about 2 cm. 


III. NEW COIL CONFIGURATION AND 
MEASUREMENTS 

The equilibration scheme we describe here uses a dif¬ 
ferent coil configuration: two coils sets of the previous 
configuration are reconnected so that the coils form an 
‘L-shaped’ coil pattern along two edges of the MSR. Four 
of these Ls are again connected in series (see Fig. |T]b). 
All connections between the coils are routed in a way 
that the fields produced by the connecting wires cancel. 
In this way, two spatial directions are equilibrated at the 
same time. The results show that the third direction can 
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FIG. 4. Simulations of the magnetization in the adjacent face 
that is not part of the equilibration loop. The geometry of 
the room and placement of the coils is the same for both 
simulations and only shown in b). In a) a current is supplied 
to the four coils at the edges, in b) all eight coils carry a 
current. The color scale is the same for both plots and in 
arbitrary units. Both cases illustrate that the currents create 
a significant magnetization in the adjacent face. 


be omitted. Only the inner layer of the MSR was mod¬ 
ified in this way. The equilibration sequence starts with 
the inner layer, followed by the unmodified outer layer 
and inner layer in sequence, using 50 s for the inner layer. 
This time can be further reduced, but below 30 s a no¬ 
table decrease in the remanent field quality is observed. 
Figure [3] shows the recorded field maps. In column a) 
the residual field for the old configuration is shown. The 
maximum value is about 0.8 nT in the lower plane and 
near the door of the MSR. In the upper planes the field 
is about 0.5 nT. The higher values close to the door are 
likely caused by stray fields of the last equilibration step. 
It was previously^ found that equilibrating the direction 
not including the door yielded the best result. However 
there stray field again magnetize the overlaps. Columns 
b) and c) show the residual field for the new configu¬ 
ration, measured twice with several days in between to 
show reproducibility. Here the residual field is slightly 
smaller, on the order of 0.5 nT or below. Also the promi¬ 
nent feature in the lower plane is not visible anymore. 
With only one equilibration cycle a re-magnetization of 
the adjacent wall is avoided. These are values that are 
achieved in typical daily operation, where non-magnetic 
furniture and tools, as well as people are inside the MSR. 

To illustrate problems with the placement of equilibra¬ 
tion coil a finite-elemente-method simulation with COM- 
SOL is shown in Fig. [I] In a) a current is applied to the 
four coils at the edges of the MSR, creating a flux loop 
around the four side walls of the MSR (around the Z di¬ 
rection). In Fig. UJd) eight coils distributed equally along 
the path of flux are supplied with half the current used 
for a). Each coil has to penetrate the top wall of the 
MSR to form a loop with shielding material enclosed. 


The simulation shows the magnetization inside the top 
wall created by the current through the coils. Since this 
face is not included in the equilibration loop, its mag¬ 
netization is not properly removed by the equilibration 
process. This dominates the error caused by a possible 
DC offset at the level of precision reached in our setup. 
Here a big number of equilibration cycles (as mentioned 
in Ref. |9) reduce this error. Additionally, when small 
volumes in the center of the MSR are considered, the er¬ 
ror has a small influence on the residual field due to the 
symmetry of the MSR. However for extended-size vol¬ 
umes these contributions become more significant. Even 
when currents are better distributed and have a smaller 
amplitude the adjacent face retains a magnetization. By 
equilibrating all faces at once this problem does not oc¬ 
cur. 


IV. CONCLUSION 

The new scheme for equilibration coils presented here 
produces static remanent fields below 0.5 nT inside a 
two-layer MSR with less contributions from magnetizable 
features at the walls. This result is on the same order of 
magnitude as that from the previous equilibration setup, 
but achieved with significantly shorter cycle times. For 
the two-layered room at TUM, the time will be reduced 
from 900 s to 150 s when both layers are changed. In an 
MSR with more layers, the reduction can be even larger. 
For the outer layers perfect equilibration is not necessary 
and the time can be shorter than the 50 s used here. If 
for example the outer three layers of a 5 layer MSR are 
equilibrated for 20 s, and the inner two for 50 s, the whole 
procedure would take only 160 s instead of 1500 s. In the 
nEDM experiment at TUM the shield has to be equili¬ 
brated for example every 10 min. With the reduced time 
for the equilibration, the duty cycle of the experiment is 
effectively increased by a factor of 2, resulting in an im¬ 
proved sensitivity of the measurement by about 40%. In 
a single layer shield, where only remanent fields of tens of 
nT are targeted, equilibration can be done within about 
10 s. Additionally, this setup only requires one set of 
coils per layer. This significantly reduces the amount of 
high current relays needed for automated switching be¬ 
tween the equilibration coils. Further, the simultaneous 
equilibration of all faces avoids accidental magnetization 
in adjacent faces due to stray fields. 
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